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The Nonperturbative Broken Phase Sphaleron Rate 

Guy D. Moore^^ 
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I present a technique for measuring the broken phase sphaleron rate nonperturbatively. There are three parts 
to the calculation: determination of the probability distribution of Chern-Simons number Acs; measurement 
of (|dAcs/dt|)jVcs=i/ 2 , the mean rate of change of Acs at the barrier; and measurement of the “dynamical 
prefactor,” the fraction of barrier crossings which result in a permanent integer change in Acs • 


1. Introduction 

The “sphaleron rate,” the spacetime diffusion 
constant for Chern-Simons number NqSj is a key 
quantity in the study of electroweak baryogene- 
sis, setting the efficiency of baryon number vio¬ 
lation. It is also a very appropriate topic at this 
conference, since it is nothing but the topological 
susceptibility of the electroweak plasma, only in 
Minkowski, time; 

r= hm ((^cs(t)-iVcs(0))b ^ 
t^oo Vt 

r can be determined reliably by perturbation 
theory in any regime where g4>{T)/T 3> a^j, 
where 4>(T) is the Higgs condensate. However, 
this is precisely the regime where it is so small 
as to be of no physical consequence. We really 
want to know it in the symmetric phase at Tc 
the phase transition temperature, to determine 
how many baryons can be produced; and in the 
broken phase at and immediately below Tc to de¬ 
termine whether they are subsequently erased or 
are preserved until the present epoch. 

Recently there has been significant progress in 
the symmetric phase rate. It has been shown 
analytically 0 and verified numerically Q that 
“hard thermal loop” radiative corrections are im¬ 
portant, and the rate is known with reasonable 
accuracy, though it has recently been demon¬ 
strated that there are logarithmic corrections [|| 
which are not well under control. 

But the real time method cannot measure the 
broken phase rate; the rate is so low that no 
events will occur in any reasonable amount of 
time. Since we cannot extend the perturbative 


calculation beyond 1 loop, we will instead go 
about a nonperturbative measurement. We can 
use the dimensional reduction approximation and 
we will relentlessly; we also take the dynamics to 
be classical, regulated to give a physically rea¬ 
sonable cutoff to the hard thermal loops. For a 
much more detailed discussion of the work pre¬ 
sented here, see [||. 

2. Defining and measuring Acs 



Figure 1. Gartoon of the free energy dependence 
on Acs 

Here is the cartoon picture of why baryon num¬ 
ber violation is slow in the broken phase. Be¬ 
cause the Higgs condensate makes being half way 
between winding numbers energetically expen¬ 
sive, there is a free energy barrier to being at 
Acs = 1/2. Our method of determining the rate 
is to find the probability distribution of Acs, and 
to multiply the chance to be in a narrow range 
about Acs = 1/2 by the mean crossing time of 
that range, set by {\dNcs/dt\). Then we need to 
determine what fraction of crossings are perma¬ 
nent. 
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In the continuum, the meaning of fVcs, modulo 
1, is the integral of FF, 

Ncs = j dr j d^xBt . (2) 

Here r parameterizes a path through the space of 
3-D configurations. The constant of integration is 
set by defining Ncs = 0 in vacuum. This is con¬ 
sistent since the above integral around any closed 
loop is the second Chern class of the loop, which is 
an integer. However, to make a nonperturbative 
study we need to find a lattice implementation for 
the above, which is touchy since no lattice defini¬ 
tion of FF is a total derivative. The best thing 
to do is to use the above definition for a unique 
specific path, the cooling path, defined by 

^ = (F“C/)(F“Fym), (3) 

dr 

with FIym the Yang-Mills part of the 3-D Hamil¬ 
tonian. This definition of Nqs is unique because 
the cooling path is; it is also minimally exposed 
to lattice problems since the cooling process very 
quickly eliminates UV fluctuations, which are re¬ 
sponsible for any poor behavior in FF. 

There are three problems we must take care 
of, though; a UV noise problem, a numerical cost 
problem, and the problem of turning a probability 
distribution for Nqs into a rate. 


3. The first two problems 


The first problem is that Nqs actually is not a 
very good measurable to be using. The problem 
is that it contains a lot of UV noise not related 
to topology. Our discussion will follow that of 
Ambjprn and Krasnitz Q. In the UV the 3-D 
classical thermal theory we are studying behaves 
at leading order like 3 copies of the abelian theory, 
in which 

^CS = I d?xe,,uF^^Al . (4) 


The mean square value at leading order in per¬ 
turbation theory is then 


(^cs) = 


g^T^V 

647r^ 


d^p 1 
(27r)3 ’ 


( 5 ) 


which is linearly UV divergent. The size of the di¬ 
vergence is extensive in volume and will be cut off 
by the lattice spacing, leading to a 1/a behavior. 

The solution is not to measure Nqs itself, but 
Ncs after some preliminary cooling has been per¬ 
formed; 

N = dr I d^xBf [F., . (6) 

To illustrate that this does the trick, we show 
Ncs and N for a period of Hamiltonian evolu¬ 
tion in the broken phase, and N for a sphaleron 
transition: 



Hamiltonian time x(g2T) 


Figure 2. Above: Ncs (oscillating like mad) and 
N (always near zero) for a Hamiltonian trajec¬ 
tory in the broken phase. Below: N during a 
sphaleron. 

This graphically illustrates that using N cuts out 
irrelevant UV noise, while still capturing interest¬ 
ing topological information. 

There is still a problem with using TV, the nu¬ 
merical cost involved in performing the cooling, 
which we must do thousands of times to get an ac¬ 
curate Monte-Carlo study of the probability dis¬ 
tribution. We finesse this problem by observ¬ 
ing that, since cooling eliminates the UV excita¬ 
tions, it quickly renders the lattice fields smooth 
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enough that we lose no information by blocking 
by a factor of 2. When we use a classically im¬ 
proved Hamiltonian and definition of FF after 
the blocking, there is almost no difference be¬ 
tween the blocked and unblocked measurement 
of -/V, provided we cool by at least r > 1.25a^ be¬ 
fore blocking. For lattices larger than about 28^ 
we can double block. 


4. Turning a probability into a rate 

We can now find the (tq dependent) probability 
distribution of N by conventional multicanonical 
means. To turn this into a rate, we first mea¬ 
sure {\dN/dt\) (sampled over iV = 1/2 configura¬ 
tions) by making a canonically distributed sample 
of fV = 1/2 configurations, generating momenta 
from the canonical ensemble, and evolving for a 
very short Hamiltonian time, measuring N before 
and after. The mean rate of = 1/2 crossings 
per 4-volume is then 

r P(|iV-0.5|<e) (|diV/dt|) ^ 

prefactor 2e V 

To measure the dynamical prefactor, which is 
the fraction of crossings which lead to permanent 
integer change in JVcs, we use the same = 1/2 
sample but continue the Hamiltonian time evo¬ 
lution, both forward and backward in time, until 
the system settles about a vacuum; then we count 
how many crossings and how many permanent JV 
changes occur. The fraction of crossings leading 
to permanent JV change is 


AT,, 
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# crossings 


1 A^fln Mnit /QN, 

0 Arg„ = W„it 


For very large Debye mass squared mf) there are 
multiple crossings, but for realistic m]y the pref¬ 
actor is about 0.4. 


5. Results 

The final results for F are compared to per¬ 
turbation theory (two loop effective potential but 
tree kinetic terms) in Figure Baryon number 
erasure after the phase transition is prevented if 
F < which occurs if a; < 0.037. 



a very large hard thermal loop strength. 



x = X/g2 


Figure 4. Nonperturbative value of F(Tc) (points 
with bars) versus perturbative estimate (curve) 
and the symmetric phase rate (curve at top), as 
a function oi x = \/g^. 
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